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We have measured the specific heat of the coupled spin-dimer antiferromagnet BaaM^Os to 
50 mK in temperature and to 29 T in the magnetic field. The experiment extends to the midpoint 
of the field region (25.9 T < H < 32.3 T) of the magnetization plateau at 1/2 of the saturation 
magnetization, and reveals the presence of three ordered phases in the field region between that of 
the magnetization plateau and the low-field spin-liquid region. The exponent of the phase boundary 
with the thermally disordered region is smaller than the theoretical value based on the Bose-Einstein 
condensation of spin triplets. At zero field and 29 T, the specific-heat data show gapped behaviors 
characteristic of spin liquids. The zero-field data indicate that the gapped triplet excitations form 
two levels whose energies differ by nearly a factor of two. At least the lower level is well localized. 
The data at 29 T reveal that the low-lying excitations at the magnetization plateau are weakly 
delocalized. 

PACS numbers: 75.30.Kz, 75.40.Cx, 75.50.Ee 



I. INTRODUCTION 

An interacting array of spin dimers, each consisting 
of two antiferromagnetically coupled spins, has a natural 
tendency to exhibit plateaus in the magnetization curve 
M(H). This is obvious if one considers such a system 
of S = 1 spins. In this case, the ground state of an 
isolated dimer is a spin singlet, the lowest excited state 
is a spin triplet, and the highest energy state is a quintet. 
An increasing magnetic field will cause the ground state 
to change from the singlet to the S z = +1 triplet, and 
then to the S z = +2 quintet, resulting in two steps in 
M(H). A weakly-coupled array of S = 1 spin dimers 
will maintain this progression of the ground state and 
will exhibit a magnetization plateau at to = 1/2, where 
to = M/M sat is the magnetization expressed as a fraction 
of the saturation value. 

For one-dimensional antiferromagnets, Oshikawa, Ya- 
manaka, and Afflecki have predicted that M(H) can have 
a plateau when the condition nS(l — m) = integer is sat- 
isfied. Here n is the period of the ground state. These 
authors have argued that the magnetization plateaus are 
generalized Haldane states with non-zero magnetization, 
implying that the transverse spin components at each 
plateau form a quantum spin liquid quite like that of the 
usual to = Haldane state. The Oshikawa- Yamanaka- 
Ameck theorem has been generalized to higher dimen- 
sions by Oshikawa 2 and refined by other authors. 3 

Plateaus at fractional values of saturation magnetiza- 
tion have been observed in Cu 2+ -based 5 = 1/2 com- 
pounds, such as the three-dimensional (3D) spin-dimer 
antiferromagnet NH4CUCI3 4 , the quasi-2D orthogonal- 
dimer antiferromagnet Sr2Cu(B03)^, and the frustrated 
quasi-2D spin system Cs2CuBr4f& NH4CUCI3 shows dis- 
tinct plateaus at to = 1/4 and 3/4, Sr2Cu(B03)2 at 



to = 1/8, 1/4, and 1/3, and Cs2CuBr4 at m = 1/3 and 
2/3. In addition, plateaus at to = 1/2 have been ob- 
served in Ni 2+ -based 5 = 1 Heisenberg antiferromagnets 
with bond alternations 

Ba3Mn 2 0s is a new 5 = 1 spin-dimer antiferromag- 
net, in which the intcr-dimcr coupling is expected on a 
structural ground to be three dimensional. The magne- 
tization of this compound has been measured at 0.65 K 
as a function of the magnetic field up to 50 T in a pulsed 
magnet»2. The data exhibit a low-field to = region and 
an m = 1/2 plateau, which occurs at around 29 T. In 
this paper, we report a determination of the magnetic 
phase diagram over the entire field range below the to = 
1/2 plateau, by means of specific heat measurements per- 
formed down to 50 mK in temperature. We also examine 
the nature of the excitations at zero field and at the to = 
1/2 plateau. 

Ba3Mn 2 0s has a trigonal structure of the space group 
R3m with the lattice parameters a — 5.711 A and 
c = 21 AAA A. The manganese ions are in the unusual 
Mn 5+ state, which makes 5=1. To date, all experiments 
including the present work have been performed on pow- 
der samples, since no single crystal has been successfully 
grown to our knowledge. The magnetization measured 
at 0.65 K is essentially zero below H c i = 9.2 T, where it 
starts to increase linearly with the magnetic fields The 
to = 1/2 plateau extends from iJ c2 = 25.9 T to H C 3 = 
32.3 T, and the saturation magnetization is reached at 
H sa t ~ 48 T. The zero-field gap A has been estimated 
from H c i to be A = g^Hci — 12.2 K. From the overall 
shape of the magnetization curve, the strengths of the 
intra-dimer and inter-dimer exchanges have been deter- 
mined to be J = 17.4 K and Ji + 2( J 2 + J 3 ) = 8.3 K, re- 
spectively, where Jo, Ji, J 2 , and J3 are the first-, second-, 
third-, and fourth-nearest-neighbor interactions, respec- 
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tively. The edges of the magnetization plateau are sharp 
in spite of the powder form of the sample, indicating an 
isotropy of the spin hamiltonian. In electron spin reso- 
nance, the absorption spectrum has a single Lorentzian 
peak with a narrow linewidth, indicating the smallness 
of any anisotropic exchange interaction and the good 
isotropy of the g factor, which is determined to be g = 
1.98i2iiS This g value has been used for the estimate of 
the gap A. The single-ion anisotropy term is estimated 
to be only D ~ 0.02 J from the linewidth, where J is the 
average of the exchange interactions^ 

II. EXPERIMENTAL METHODS 

The powder sample of BaaM^Os used in the present 
work was synthesized by the method described in Ref.lTll 
with an extended sintering time of 150 h to improve the 
sample quality. 9 Specific heat measurements to 18 T were 
performed in a superconducting magnet at temperatures 
down to 50 mK using a dilution refrigerator. The re- 
laxation calorimeter for this setup has been described in 
Ref. 0. For measurements from 19.5 T to 29 T, an- 
other relaxation calorimeter with a built-in 3 He refriger- 
ator was used in a resistive Bitter magnet. 

Powder samples present a challenge to relaxation 
calorimetry at temperatures typically below 1 K be- 
cause of the difficulty in obtaining a good thermal con- 
tact between the powder grains and the calorimeter. To 
solve this problem, we mixed the sample with a small 
amount of silver painlji^ and sandwiched it between two 
0.13 mm-thick pieces of sapphire to form a thin, uniform 
layer. The bottom piece of sapphire was attached with 
the Wakefield compound^ to the calorimeter platform. 
Three samples weighing 0.15 mg, 1.08 mg, and 7.62 mg 
were used to cover different temperature and field re- 
gions. 

III. EXPERIMENTAL RESULTS AND 
DISCUSSIONS 
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FIG. 1: Temperature dependence of the specific heat of 
Ba3Mn2 08 in constant magnetic fields between 9 T and 18 T. 
The nuclear contribution of 55 Mn is visible at temperatures 
below 0.2 K. 
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FIG. 2: Specific heat of BasM^Og in constant magnetic 
fields ranging from 19.5 T to 26 T. The lines are guides to the 
eye. 



A. Magnetic phase diagram 

The temperature dependence of the specific heat is 
shown in Fig. 1 for fields between 9 T and 18 T. In fields 
higher than 13 T, a well defined peak appears in the spe- 
cific heat below 1 K, signalling a phase transition. The 
sharpness of the peak indicates a high quality of the sam- 
ple despite its powder form. An additional peak which 
indicates a second transition appears at 12 T. In fields 
lower than this value, the specific heat anomaly becomes 
weak, and we can identify only one peak with confidence. 
No peak is found at 9 T. 

The specific heat in magnetic fields higher than 18 T, 
measured in the resistive Bitter magnet, is shown in 
Fig. 2. The peak in the specific heat is sharp up to 22 T. 
The 24 T data are very similar to the 12 T data, except 



that there are two shoulders instead of two peaks. At 
25 T, the anomaly becomes small, and only one broad 
peak can be identified. 

To further explore the field region below 13 T, we have 
measured the specific heat as a function of the magnetic 
field. In these measurements, a constant electric current 
was applied to the heater of the thermal reservoir, allow- 
ing the temperature to rise monotonically with increasing 
field as dictated by the magnetoresistance of the heater. 
As seen in Fig. 3, the phase transitions are observed as a 
peak and shoulder. 

The magnetic phase diagram determined from the po- 
sitions of the specific-heat peaks and shoulders in Figs. 1 
through 3 is given in Fig. 4. As seen in the figure, the 
transition temperatures obtained in the magnetic-field 
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FIG. 3: Magnetic-field dependence of the specific heat of 
Ba3Mn2 08- The transition points are indicated by arrows. 
The lines are guides to the eye. 



scans shown in Fig. 3 agree with those obtained in the 
temperature scans shown in Fig. 1. The phase diagram 
contains three ordered phases in the investigated field re- 
gion. In addition to the H c \ and H C 2 observed in the 
magnetization dataS, there are two critical fields at the 
intermediate fields of about 11 T and 24.5 T. Figure 4 
is somewhat reminiscent of the magnetic phase diagram 
of CsFeBr3(A£i which is an S — 1 antiferromagnet on a 
triangular lattice with a singlet ground state. In this ma- 
terial, H c of the field-induced ordering at zero tempera- 
ture splits into two branches as the temperature is raised, 
with the separation between the two branches widening 
with increasing temperature^ In contrast, Fig. 4 sug- 
gests two branches of either H c i or H C 2 that are separate 
at zero temperature merging into one as the temperature 
is raised. 

The Bose-Einstein condensation of triplets is believed 
to be a valid description of field-induced long-range or- 
dering in all gapped antiferromagnets^ii^ According to 
the Bose-Einstein condensation theory, the phase bound- 
ary in the temperature vs field plot obeys a power law 
T c oc (H — H c ) a . The field-induced ordering in the 
5=1/2 spin-dimer material TICUCI3 has been inter- 
preted in terms of the Bose-Einstein condensation^ 
However, the exponent a of 0.50 determined for this ma- 
terial is smaller than the theoretical value of 2/3 obtained 
with a Hartree-Fock approximation. 

We have obtained H cl = 9.04 ±0.15 T and a = 
0.39 ± 0.06 by fitting the data of the phase boundary 
between the disordered phase and the ordered phase be- 
low 0.58 K. The zero-temperature critical field H c i is in 
good agreement with 9.2 T obtained by the magnetiza- 
tion curve at 0.65 K. Note that the critical fields deter- 
mined by a low-temperature magnetization curve M{H) 
are generally closer to the zero-field values than the ac- 
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FIG. 4: Magnetic phase diagram of BaaM^Os- The circles 
are the positions of the peaks and shoulders in the specific 
heat at constant magnetic fields. The diamonds are the posi- 
tions of the peaks and shoulders in the specific heat measured 
as a function of the magnetic field, while the current for the 
reservoir heater was held constant. The lines are guides to 
the eye. 



tual transition fields at the temperature of the measure- 
ment. The value for the exponent a is smaller than 2/3 
predicted by the theory. However, it agrees within the 
experimental uncertainty with the quantum Monte-Carlo 
result a = 0.37 ± 0.03 for S — 1/2 spin dimers with a 
weak three-dimensional inter-dimer coupling. 19 It is also 
close to the values for the S = 1 spin-chain materials 
NTENP and NDMAP for the field applied parallel to 
the chain direction, a = 0.334 ± 0.004 and 0.34 ± 0.02 
for NTENBSa and NDMAPjSi respectively. Recently, 
careful analyses near the zero-temperature critical field 
H c i have shown that the exponent a approaches 2/3, as 
the temperature region of the fits is narrowed, for both 
the quantum Monte-Carlo data2£ and the experimen- 
tal dataS for TICUCI3. A similar analysis has yielded 
a = 0.63 ± 0.03 for BaCuSi 2 6 , an S = 1/2 spin- 
dimer compound with predominantly two-dimensional 
inter-dimer exchanges^ It remains to be seen, however, 
whether the Bose-Einstein condensation theory can ex- 
plain the field-induced magnetic ordering in all gapped 
antiferromagnets. 



B. Gapped excitations at zero field and at the 
magnetization plateau 

In the low-field m — region, as well as at the to = 1/2 
magnetization plateau, the specific heat is expected to ex- 
hibit a gapped behavior. The data at zero field and 29 T, 
the midpoint of the plateau, are shown in Fig. 5. At zero 
field, the low-temperature behavior of the data indicates 
the presence of a gap. However, the low-temperature spe- 
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FIG. 5: Temperature dependence of the specific heat of 
Ba 3 Mn 2 8 at T and 29 T. The thick solid lines are the 
low-temperature fits to the data. The dashed line indi- 
cates the calculated zero-field specific heat due to localized, 
gapped excitations, whose parameters are taken from the low- 
temperature fit. The thin solid line is the specific heat due 
to two kinds of localized excitations with two gap energies, 
as described in the text. The dotted line for the 29 T data 
indicates the specific heat of mutually isolated dimers with 
an excitation gap of A = 4.3 K. In the inset, ln(CT 2 ) of the 
zero-field data have been plotted against 1/T to reveal the 
localized-excitation behavior, which is given by the straight 
line. The deviation of the data from the line at temperatures 
below 1.4 K is due to the phonon contribution and an uncer- 
tainty in the subtraction of the background heat capacity of 
the silver paint. 



cific heats of delocalized gapped excitations have wrong 
temperature dependences — see Eq. 4 for an example — 



which do not agree with the data. The best fit of the data 
below 3.5 K is obtained by using the low-temperature for- 
mula for localized gapped excitations, 



C = ni?(A/T) 2 e- A / T , 



(1) 



where n is the number of excited states per spin dimer, 
and R is the gas constant. As shown in the inset of 
Fig. 5, this formula gives the best fit with an energy gap 
ofA = 12.5±0.2K, which is in excellent agreement with 
12.2 K estimated from H&. In BasM^Os, the average 
inter-dimer exchange J\ + 2( J2 + J3) = 8.3 K is substan- 
tial in comparison with the intra-dimer exchange Jo = 
17.4 K, as we have described. It is surprising that, de- 
spite such strong inter-dimer exchanges, the triplet exci- 
tations remain localized instead of becoming highly dis- 
persive. The crystal structure^ of BasM^Os suggests 
that the localization arises from a geometric frustration. 

Another surprising result of the low-temperature fit is 
that n is only 1.44 ± 0.04, which is very close to 1/2 of 
n = 3 expected for the triplets. The dashed line in the 
main panel of the figure shows the contribution of local 
excitations for the entire temperature range of the data, 
according to the complete expression 



C = hR 



(A/T) 



2 p -A/T 



(l + ne- A / T ) 2 ' 



(2) 



for h — 1.44 and A = 12.5 K found by the low- 
temperature fit. Above 3.5 K, the specific heat is sub- 
stantially larger than can be accounted for by these ex- 
citations, strongly suggesting that the excess arises from 
the "missing" triplets whose gap is considerably larger 
than 12.5 K. Indeed, as shown by the thin solid line, the 
data over the entire temperature range can be well de- 
scribed by 1/2 of the triplets at A = 12.5 ± 0.2 K and 
the other 1/2 at A' = 24.5 ± 0.5 K: 



J 



. (A/T) 2 e- A / T + (A'/r) 2 e - A '/ T + fi(A'/T - A/r) 2 e -( A+A '>/ T 
[l + n(e- A / T + e- A '/T)]2 



(3) 



r 



where n — 1.5. We find A' = 2 A within the combined 
uncertainties of the fitting parameters. Although Eq. 3 
assumes that the excitations at A' are localized like those 
at A, the data do not tell whether this is the case, since 
the excitations at A dominates the low-temperature be- 
havior which distinguishes localized excitations from dis- 
persive excitations. 

The field of 29 T is located at the midpoint of the 
m = 1/2 magnetization plateau within the uncertainty 
of the magnetization data. As expected, no sign of or- 
dering is observed in the specific-heat data at this field, 
indicating that the plateau region is a spin-liquid phase, 
although there is a substantial scatter at high temper- 



atures. At the midpoint of an m = 1/2 magnetization 
plateau, mutually isolated S = 1 dimers will have the 
lowest excited level containing two states, a singlet and 
an S z = 2 quintet. The gap for the level of the col- 
lective excitations analogous to such single-dimer states 
can be estimated from the width of the magnetization 
plateau to be A = g^,Q(H c3 — H c2 )/2 — 4.3 K. The spe- 
cific heat due to localized gapped excitations with h = 
2 and A = 4.3 K, shown by the dotted line, completely 
fails to describe the data at 29 T. The comparison of the 
line with the data in fact indicates that the energy of the 
excitations has a wide distribution, suggesting that their 
dispersion cannot be neglected in comparison with the 
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gap size. 

By following the method described by Troyer et alw2> 
for one-dimensional excitations, one can derive the low- 
temperature expression 



C = n m R 



A 2 



(4 7 ra) 3 / 2 T 1 /2 



-A/T 



(4) 



for the specific heat due to three-dimensional excitations 
whose gapped dispersion is e = A + a(fc — ko) 2 . Here n m is 
the number of the gapped modes. This expression gives 
the best fit to the data below 1.5 K as shown in Fig. 5, 
with A = 1.9 ±0.1 K, which is considerably smaller than 
4.3 K anticipated by the simple argument given above. 
The strength of the dispersion can be estimated from 
the fitting parameter a/rim 3 — 0.20 ± 0.01 K. Regard- 
less of the number of modes n m , which we do not know 
at the moment, this is a rather small value, indicating 
that the low- lying excitations at the m = 1/2 magnetiza- 
tion plateau are nearly localized, similar to the low-lying 
excitations at zero field. The measured specific heat con- 
tinues to rise above 1.5 K, where the low-temperature ap- 
proximation is invalid, revealing the presence of at least 
another low-lying mode whose energy is on the order of 
10 K. This observation, together with the experimentally 
obtained value for A being considerably smaller than the 
value estimated from the plateau width, suggests that the 
levels of the singlet and the S z — 2 quintet do not cross 
each other at the midpoint of the magnetization plateau. 



IV. CONCLUSIONS 

In this work, we have determined the magnetic phase 
diagram of BaaM^Os by using specific- heat measure- 
ments to 50 mK in temperature and to 29 T in the mag- 
netic field. The phase diagram reveals the existence of 
three ordered phases between the low-field spin-liquid re- 



gion and the m = 1/2 magnetization plateau. We have 
found that the exponent for the transition temperature 
of the field-induced ordering in the investigated tempera- 
ture region is clearly less than 2 /3 expected by the Bose- 
Einstein condensation theory. 

In the spin-liquid phase at zero field, the specific heat 
over the entire temperature range is well accounted for by 
two excitation levels, each probably containing the same 
number of states. The gap A for the lower level agrees 
well with the value estimated from H&, and the excita- 
tions at this level are well localized. The existence of the 
higher level at about 2A requires direct confirmation by 
ESR or inelastic neutron scattering. 

In the spin-liquid phase at the m = 1/2 magnetiza- 
tion plateau, the low-lying excitations are delocalized, 
but the dispersion is weak. Contrary to a naive expec- 
tation, the two collective modes arising from the singlet 
and the S z = 2 quintet of a single dimer do not cross 
each other at the midpoint of the plateau. We hope this 
finding will stimulate further theoretical work on three- 
dimensionally coupled spin dimers. In addition, it will 
be useful to perform an ESR experiment to map the ex- 
citation energies in the magnetization plateau region as 
a function of field. 
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